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Abstract. The combination of stellar population synthesis models and photoionization models 
allows a better understanding of the spectral features of HII regions and HII galaxies. In this work 
we show that sampling effects in the initial mass function (IMF) are very important in the low 
cluster mass case. To this aim, we compute photoionization models ionized by realistic clusters 
- - - , made up of various combinations of individual stars and clusters made up with a synthesis 
■ model. We discuss the differences in the position on diagnostic diagrams and their implications. 
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1. Introduction 



I ■ Clusters modeled by stellar population synthesis models have been used as photoioniza- 
[ tion sources in several works (|Cerviiio fc Mas-Hesse 19941ILuridiana et 'al. 1999 , Stasinska fc Izotov 20031 
T-H . among others). Most of the current models assume that the Initial Mass Function (IMF) 
^ ' is completely sampled, but when the number of stars is small this assumption is in- 
■ correct. Therefore standard synthesis models cannot be applied to low mass clusters. 
\ Cervino et al. (2003) Cervino & Luridiana (2004 2006) showed that this question can 
' be treated from a probabilistic point of view in which the IMF is interpreted as the prob- 
ability distribution function of the stellar mass rather than its deterministic expression. 
C ■ Implementing the formalism developed by the mentioned authors, the sampling problem 
can be studied in a better way. The aim of this work is to show the importance of sam- 
pling effects in the IMF for low mass clusters, the most abundant ones, and to emphasize 
the need of a probabilistic approach to this problem. 

o ■ 

c/3 , 2. Method and results 

\ Using Emission Line Diagnostic Diagrams (ELDDs) we compare the emission-line spec- 
^ ■ tra obtained with a synthesis model (Cervino et al. 20021) with others obtained with com- 
k>( , binations of individual stars of 20, 25, 50 and 100 Mq. In each of these combinations 
5_j ■ the stars have the same mass so that sampling effects are important. With these models 
we can explore a three parameter space: mass, number of ionizing photons (Q(H")) and 
number of stars. We consider three cases: fixed mass equal to 100 M0; fixed Q(H°) equal 
to the one of a 100 M0 star; and fixed number of stars equal to N=:l. All the models 
have been computed for ZAMS. 

Using CLOUDY (jFerland et al. 1998P we compute the emission of model nebulae with 
Orion abundances ionized by the clusters described above and calculate their represen- 
tative points in ELDDs. 

Figure 1 (left) show the differences in a selected ELDD between the predictions of the 
synthesis model cases and those of the rest of the cases. The synthesis model cases fall 
either far from any of the stellar cases or very close to the 50 Mq stellar case. Moreover, 
the points corresponding to the synthesis model cases closely follow the tendency of the 
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Figure 1. Left: [O III] A5007/H/3 vs. [O I] A6300/Hq intensity ratios, dashed lines join cases 
with stars of the same mass. Right: ionizing continua for fixed Q(H'') cases. 

50 M0 cases. This is expected because the ionizing continuum for the synthesis model 
case is very similar in shape to the one 50 Mq star (Fig. 1, right). In all the cases the same 
abundances are used, but there is a considerable dispersion in the position on ELDDs 
due to the differences in the ionizing continua of the stars used. So, if a comparison with 
theoretical grids of HII regions based on usual synthesis models were done, it could be 
wrongly deduced that our theoretical cases span a wide range in metallicity and/or age. 

3. Conclusions 

The dispersion shown by empirical ELDDs can be attributed not only a spread in 
metallicity and age but also to IMF sampling effects, which introduce an additional 
degeneration. The smaller the cluster, the more severe the effect. If we can infer from 
observations the value of Q(H'') or the mass of the underlying ionizing cluster it's impos- 
sible to accurately determine, by means of a standard synthesis model, the corresponding 
point in none of the ELDDs given the separation between the models. Therefore stan- 
dard synthesis models should not be applied to low mass clusters because they would 
lead to misinterpretations of observational or theoretical results. In order to comple- 
ment this work, the relative probabilities of the different cases will be computed in the 
next future using the probabilistic formulation followed by Cervifio & Luridiana (|2006)[ 
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